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ABSTRACT 

The observed distribution of galaxies has local transverse isotropy around the line-of- 
sight (LOS) with respect to the observer. The difference in the statistical clustering 
signal along and across the line-of-sight encodes important information about the ge¬ 
ometry of the Universe, its expansion rate and the rate of growth of structure within 
it. Because the LOS varies across a survey, the standard Fast Fourier Transform (FFT) 
based methods of measuring the Anisotropic Power-Spectrum (APS) cannot be used 
for surveys with wide observational footprint, other than to measure the monopole 
moment. We derive a simple analytic formula to quantify the bias for higher-order 
Legendre moments and we demonstrate that it is scale independent for a simple sur¬ 
vey model, and depends only on the observed area. We derive a similar numerical 
correction formula for recently proposed alternative estimators of the APS that are 
based on summing over galaxies rather than using an FFT, and can therefore in¬ 
corporate a varying LOS. We demonstrate that their bias depends on scale but not 
on the observed area. For a quadrupole the bias is always less than 1 per cent for 
k > O.OlhMpc”^ at z > 0.32. For a hexadecapole the bias is below 5 per cent for 
k > 0.05hMpc”^ at z > 0.32. 

Key words: methods: data analysis — methods: numerical — galaxies: statistics — 
dark energy — distance scale — large-scale structure of Universe. 


1 INTRODUCTION 


The three-dimensional Power Spectrum (PS) of galaxies is 
one of the most important measurements that can be made 
from galaxy surveys. The Baryon Acoustic Oscillation fea¬ 
ture in the PS can be used to obtain sub percent constraints 
on the expansion history of the Universe; The Redshift- 
Space Distortions (RSD; [Kaiser I pST] ) allow precise mea¬ 
surements of the growth rate of structure; And the Alcock- 
Paczynski (AP; Alcock fc Paczynski||1979 1 effect constraints 
very tightly the geometry of the Universe (for the most re¬ 
cent measurements see, e.g. [Anderson et al.|2014[ jSamushlal 
et al. 12014 1. 


Both RSD and AP are imprinted into galaxy distribu¬ 
tion as a signature along the line-of-sight (LOS) from the 
observer. To extract these signals in an unbiased way it is 
important that we analyse the data using the correct LOS 
that varies from a galaxy-pair to a galaxy-pair. It is geomet¬ 
rically impossible to define a Cartesian coordinate grid in 
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such a way that the z-axis is everywhere aligned with the 
LOS direction. Thus, the APS cannot be measured in Carte¬ 
sian coordinates with one of the directions in dual Fourier 
space identified with the LOS. 


Because the PS is a 2-point statistic, it relies on the 
properties of pairs of overdensities, although estimation 
methods can instead be based on pairs of galaxies. The vary¬ 
ing LOS means that, for pairs of galaxies separated by wide- 
angles, the RSD for the galaxies in a pair will not be parallel. 
The resulting clustering signal including these wide-angle 


(WA) effects can be accurately modelled (Szalay, Matsub- 
ara, & Tandy 1998 Szapudi 20041, but the difference beyond 


assuming a single LOS to the mid-point between the pair 
of galaxies is small (Samushia, Percival, Sz Raccanelli[[2012 


Beutler et al. 2012 Yoo & Seljak 20151. In this paper we 


will concentrate on quantifying the effect of different meth¬ 
ods to allow for the varying LOS between different pairs of 
galaxies. We will assume that the WA effects are small. 


For distant surveys covering a small spatial region, the 
LOS will not vary significantly across the survey. In order to 
see where the approximation of a single-LOS breaks down, 
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suppose that we consider making a general FT of the over¬ 
density field, 

5{k) = J d^r (5(r)e~‘^’’, (1) 

using an FFT algorithm, with z-axis pointed towards the 
middle of the survey. The APS could then be computed by 
averaging S{k) in wavenumber bins, 

P{k) = ^(5ik)S*ik)), ( 2 ) 


where V is the volume of the survey. This method has been 
successfully used in the past to measure the monopole mo¬ 
ment of the power spectrum, where the LOS direction is 
irrelevant, and we will now contrast it with recent methods 
for measuring higher-order moments. We will refer to this 
method as a single-LOS method. 

What we actually want is a method to compute the 
APS, or moments of it, as if there were no LOS variations, 
i.e. automatically correct for LOS variations in the compu¬ 
tation of the APS, allowing the fast comparison with mod¬ 
els and the retention of all information. The most natural 
way of doing this, allowing for a radially-orientated LOS is 
to perform an FT in spherical coordinates with the origin 
at the observers position. Even though approach has been 


applied to previous galaxy surveys (see, e.g. Tadros et al. 


1999 Percival et al.|[^004 Leistedt et al.|[^012 l, FFTs have 

an advantage of being faster and more efficient in terms of 
computing both the PS estimators and the subsequent like¬ 
lihood analysis of the measurements. 

In principle, it is also possible to measure the corrected 
APS on a Cartesian grid by replacing the FFT with a sum 


over galaxy-pairs (Yamamoto et al. 20061. This allows a 


LOS to be defined for each galaxy pair in the sample (ef¬ 
fectively having a separate coordinate grid, or z-axis for 
each galaxy pair). It is however impractical for modern large 
galaxy surveys, as it would require prohibitively large com¬ 
puting times. We will refer to this method as a pairwise-LOS 


method (see Sec. 2.11. As discussed above, in this paper we 
will ignore WA effects, and thus consider that the pairwise- 
LOS method gives an exact result. 

A practical approximation for the pairwise-LOS method 
is to define a LOS for a chosen galaxy in each pair. This al¬ 
lows the estimation method to be reduced from a sum over 
pairs to a sum over galaxies, which is computationally faster 
I Blake et al. pm| . This approximation will break down for 
galaxy pairs with very large angular separation (it is ef¬ 
fectively another WA effect) but will become increasingly 
accurate for smaller scale measurements. The algorithm can 


be reduced to series of FFTs (Bianchi et al. 2015 Scocci- 
|marro | [M5| and is significantly faster than the pairwise- 
LOS algorithms, which makes it feasible for the analysis of 
galaxy surveys. This method has been applied to WiggleZ 
data ( [Blake et al.|20lf| and BOSS data ( jBeutler et al.|20T^ 
and the results suggest that there isn’t an appreciable bias 
with respect to the pairwise-LOS method. We will refer to 
this method as a moving-LOS method (see Sec. |2.3| ). 

In this paper we aim to quantify the biases in the APS 
induced by single-LOS and moving-LOS methods with re¬ 
spect to the pairwise-LOS method. To make the discussion 
clearer we make certain simplifying assumptions: We ignore 
the effects of mask and selection function (window effects). 


Correcting for the mask and selection effects is not a trivial 
task, but this problem is almost independent from the issue 
that we want to address in this work, so we will assume that 
the window effects have been properly dealt with to required 
accuracy. We also ignore the discrete nature of galaxy survey 
data and will write all equations as integrals over overden¬ 
sity field rather than sums over galaxies. These assumptions 
help to keep the discussion clearer and the equations com¬ 
pact and don’t affect any of our main conclusions. The ap¬ 
proximation is further justified by the fact that for large 
separations (where the geometric bias is larger) the effects 
of discreteness of the galaxy field are negligible. 

We denote vectors by bold symbols (r), unit vectors by 
bold symbols with a hat (r = r/|r|), and the modulus of a 
vector with italic symbols (r = rr/\r\^). A scalar product is 
assumed between two sequential (unit) vectors. 


2 ANISOTROPIC POWER-SPECTRUM 

We start with the basic premise of a correlated galaxy over¬ 
density field, with the correlation function (CF) between two 
galaxies at positions ri and r 2 , ^{ri, r 2 ) = {S{ri)5{r2)}. For 
convenience, we define two vectors 

r- = r 2 -ri, (3) 

r+ = (r2-|-ri)/2, (4) 

where r- connects the two galaxies and r+ is the vector from 
the observer to to their midpoint, which we will identify 
with the LOS of the galaxy pair. Because of local transverse 
isotropy around the line-of-sight, and our assumption of no 
WA effects, the CF will only be a function of the distance 
between the galaxies and the angle they make with respect 
to the LOS. ^(ri,r 2 ) = ^(r-,r-r+). 

The angular dependence of the CF is usually expanded 
into Legendre polynomials 

^{r-,r-r^) = (5) 

e 


with most of the useful information in first three even mul¬ 
tipoles (Taruya, Saito, fc NishimichipOlT Kazin, Sanchez, 


jfc Blanton|2012yr 

The APS is defined as a FT of CF and can also be 
decomposed into Legendre polynomials with respect to LOS 


P(fe) = / d^r-^(r-, r-T%)e ^ P(!(fc)Tr(fcr+). (6) 


This is a standard definition and theoretical predictions of 


APS are computed for this quantity (see e.g., Reimberg, 
Bernardeau, & Pitrou 2015 and references therein). The 
PS multipoles are related to the CF multipoles by 


Pe{k) = ^rri^y dr- ^e(r-)je(kr-)r^. 


( 7 ) 


2.1 Pairwise-LOS Method 

In the pairwise-LOS method one would correct for the vary¬ 
ing LOS, by computing the integral over the overdensity 
field simultaneously assigning correct LOS direction to all 
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(a) (b) 

Figure 1. (a) Angles between the pair separation vector and the LOS assumed in the pairwise-LOS method (red) and the single-LOS 
method (green), (b) Angles between the pair separation vector and the LOS assumed in the pairwise-LOS method (red) and the moving- 
LOS (green) methods, for two example pairs of galaxies AA’ and BB’. The moving-LOS method fails (i.e. red and green angles are 
significantly different) for pairs with large separations, while the single-LOS method fails for pairs whose true LOS is different from the 
fixed LOS assumed in the single-LOS method. The dashed arrows on panel (b) are parallel to lines OB and OA. 


galaxy pairs. In other words, one would compute the multi¬ 
dimensional integral 

Jdkd^rid^r2 

6{ri)6{r2)e~''^'^^ Ce{k7\) ■ ( 8 ) 

The expectation value of this integral is 

PAk) = f 

d-jiV J 

^(r-, r-f+)e~‘^’'’£^(fcr%) (9) 

and the APS computed in such way would coincide with the 
definitions of Eq. § and 0 and would therefore be exact. 
Computing this multi-dimensional integral over large galaxy 
sample however demands prohibitively large CPU time and 
is not currently viable. 


2.2 Single-LOS Method 


In the single-LOS method the z-axis is pointed towards the 
middle of the survey footprint and it’s assumed that r + ~ £ 
within the survey volume. This approximation allows us to 
rewrite Eq. ^ as 


2/4-1 f ^ ^ 

Pdk) = J dkCeikz) 



( 10 ) 


The expectation value of the integral is 


Pdk) = JdkPik)£i{kz), ( 11 ) 


where 

P{k) = ^ JJd^r-d^r+Gr-,P-P*)e~^'^'^~, (12) 

as before. After integrating over r- and k this expression 
reduces to 

Pi{k)^P^Gk)^ jdG.Ceizr,), (13) 

where P* is the true PS multipole defined by Eqs. 
and Q. 

For / = 0 the single-LOS APS reduces to the true APS 
as £-o{x) = 1 - The monopole of the APS is unbiased even 
when measured with the single-LOS method. For higher or¬ 
der multipoles the bias is always present, but could be small 
in some limits. For example, if the survey area is small and 
the z-axis is pointed towards the center of the survey, LOS 
directions of all pairs will be very close to z. In this case, 
r+ ~ z and £r(zr+) ~ 1, resulting in small bias. 


2.3 Moving-LOS Method 

The moving-LOS method is an approximation of Eq. 

The overdensity field is transformed as 

5i{k) = jd?r 5{r)Ci{kr)e~^'^'" (14) 

and the APS multipoles are computed as 

Pt{k) = Jdk5oik)5i{k). (15) 
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The APS computed in using this expression is equiva¬ 
lent to 

Pi(k) = J d^rid^r2dfc 

5(r 1 ) 5(r2 )e-‘'“"i Ce (kf]_) (16) 

and the expectation value of that expression is 

Pe{k) = } f d^r-d^r+dfc 

a-rV j 

^(r-, f-f+)e“‘^’’'£r(fefi). (17) 

Using properties of spherical harmonics this can be further 
reduced to 

Pr(fc)= f dk'Y,PHk')Wu>{k,k'), (18) 

■’ i' 

where 

2 t '2 r 

Wt^t{k, k') = —— J dr-r^je{kr-)je,{k'r-)Pu,{r-) (19) 
and 

Yi>m' {r-)Yi,^, {r-)Yem (fi) 

mm' 

(See App. I^for details). 

This expression will reduce to the true APS only if 
Wu'{k, k') = Seei5{k — k'). In the limit of very small separa¬ 
tions r+ ~ ri, the P tends to 5u' by the virtue of the closure 
relation for the spherical harmonics, making Wu'{k, k') con¬ 
verge to Su'&(k — k') in the meanQ Properties of spherical 
harmonics also enforce the condition Wor' oc Sot' for all k 
and k', which means that, as with the single-LOS method, 
there is no bias for £ = 0. Unlike the single-LOS method, 
the bias can’t be expressed as a simple scale-independent 
scaling of the APS. 


3 APS BIAS AS A FUNCTION OF SKY AREA, 
REDSHIFT, AND SCALE 

In this section we will quantify the biases in the single-LOS 
and moving-LOS methods, using the pairwise-LOS method 
as the reference. To compute the biases we need to specify 
the geometry of the observed volume as the biases will de¬ 
pend on the distribution of pair separations. For simplicity, 
we will assume that the observed patch of the sky is circu¬ 
larly symmetric, the z-axis is pointing to the center of the 
observed area (this choice results in the least bias for the 
single-LOS method with the assumed LOS along this direc¬ 
tion), and the mask and selection functions are uniform. We 
will also assume that the width of redshift bin is small com¬ 
pared to the distance of the sample from the observer. This 
simple model clearly lacks the detailed window of an actual 
survey, but the angular distribution of pair separations will 


^ In this limit the moving-LOS method basically reduces to the 
pairwise-LOS, since the pairwise-LOS method will have the same 
expansion as Eq. IpOl but with ri replaced by r+. 



Figure 2. Fractional error in the measured APS as a function of 
observed area for the single-LOS method. 


be roughly correct, and the radial thinness is a conserva¬ 
tive choice as it forces pair separations of the same physical 
separation to wider angular separations. 

For this geometry, using a single-LOS analysis the bias. 


given by Eq. (131, reduces to 


Ptik) = P^k) 


f 

Jo 


di? sin(T?)£^[cos(i?)] 


/ 


( 21 ) 


di9 sin(d) 


The bias in the APS is independent of the wavenumber and 
redshift and only depends on the angular extent of the ob¬ 
served area. We plot this bias for 1 = 2 and I = A m 
Fig.m We see a gradual increase in the biases with area 
for small surveys, with biases that are already larger than 
1 per cent when the footprint is only 1000 deg^. For hemi¬ 
spherical or full-sky surveys with a single-LOS, both the 
quadrupole (I = 2) and the hexadecapole {I = A) are zero 
(the fractional error is 1): here any increase in the cluster¬ 
ing strength caused by RSD along the single-LOS is matched 
by an increase perpendicular to the single-LOS around the 
edges of each hemisphere. 

For the moving-LOS method, the APS bias given by 


Eq. (181 cannot be expressed as a simple ratio of true and 


measured power-spectra. For the simple - “thin spherical 
cap” - geometry we have assumed, we can use the proper¬ 
ties of spherical harmonics to reduce the five dimensional 


integral in Eq. (201 to a one dimensional integral 


Pw iv) =1 


.t-t' 2£-bl 


i: 


dp 


( 22 ) 




1 - 7?p/2 


yi - ryp -b r;2/4 


(see, App. 0 where P denote associated Legendre poly¬ 
nomials and rj = r-jr^. Examining Eq. ( |22[ ) we see that 
the expression tends to Su' when r] tends to zero, suggest¬ 
ing that, for a fixed scale, the approximation works better 
the further away that galaxies are from the observer, as ex¬ 


pected. Eq. (221 also shows that the bias depends only on the 
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wavenumber (h/Mpc) 


Figure 3. Fractional error in the measured APS as a function 
of observed area for the moving-LOS method for z = 0.32 and 
0 = 0.35. 


wavenumber and doesn’t depend on the observed sky cover¬ 
age: given that the bias is only related to how each pair of 
galaxies is treated rather than the distribution of pairs, this 
is also expected for scales unaffected by the window. 

Fig. [3] shows the fractional bias in the APS at z = 0.32 
for /3 = 0.35[^The bias in the quadrupole lies below 1 per 
cent at wavenumbers above 0.01 h/Mpc. The bias in the 
hexadecapole is larger and reaches a sub per cent level only 
for wavenumbers above 0.1 h/Mpc. At higher redshifts the 
biases are reduced even further. 


SUMMARY AND COMPARISON WITH 
PREVIOUS WORK 


In this work we quantify biases on the APS measured by var¬ 
ious computational methods, focussing on methods correct¬ 
ing for the varying LOS. We have ignored wide-angle effects 
that arise because the peculiar velocity shifts in galaxies in 
a pair are not parallel. At 200/i'^ Mpc this effect is less than 
a per cent already at 2 = 0.2, and decreases for the smaller 
scales usually of interest in analyses (see e.g.|Samushia, Per 


cival, &: Raccanelli||2012] |Raccanelli, Samushia, fc Percival 


2010). Instead, the approximations discussed in this work 


are related to the way of measuring the APS, and we have 
compared two methods - one known to be wrong (single- 
LOS), and one where small biases have previously been as¬ 
sumed ( [Blake et al.||20lT| [Beutler et al.|[^014| ) but that we 
now quantify (moving-LOS). 

To provide a baseline, we have considered ignoring the 
variation in the LOS and assuming a single-LOS: here we 
have shown that we are left with a significant bias in mea¬ 
surements other than the monopole moment of the power 
spectrum, that is independent of scale and only depends on 
the area of the angular footprint. This has been know for 


^ 2 = 0.32 is an effective redshift of BOSS CMASS sample. /3 
parameter describes the amount of anisotropy in APS (see, e.g. 


Hamilton|1998 for a proper definition). 


a long time, and consequently this method has not been 
widely used to measure the APS, other than the monopole. 
We have shown through an analytic formula for the in¬ 
duced bias of APS multipoles (Eq. 13), that the bias on the 
quadrupole is large even for surveys as small as 1000 deg^. 
Our results can be used to correct this effect to leading or¬ 
der for analysis methods in which the observed area is sub¬ 
divided into smaller regions to make the bias smaller (see 
e.g., [Cole, Fisher, fc Weinberg||1904l [Hemantha, Wang, fcj 
Chuang|2014). 


We have also considered a measurement method pro¬ 
posed more recently, that performs the transform summing 
over galaxies in a way that allows varying LOS to be factored 
into the measurement. Ideally we would want to sum over 
pairs rather than galaxies, but this is impractical, and this 
revision leads to a small bias in the anisotropic measure¬ 
ments. We have derived an analytic formula for this bais 
that, in the limit of narrow redshift bin, can be reduced 
to a one dimensional integral over the true power-spectrum 
(Eq. |18| ). This shows that the bias is a function of scale 
(but not the area) and only depends on the ratio of pair- 
separation to the distance of pair from the observer r-/r+. 
While small, this correction can easily be calculated and 
measurements corrected for this effect. 

The new method calculating the APS using a sum over 
galaxies (moving-LOS) remains accurate for the quadrupole 
APS at scales above k — Mpc even at small red- 

shifts. The bias in the hexadecapole moment is larger, but 
also decays at high redshifts. Our correction formulas in 
Eq. (13) and (181 can be used to correct the biases to leading 
order. 

Our analysis of the significance of making various LOS 
approximations makes a number of simplifying assumptions 
about the survey for the sake of analytical clarity. We as¬ 
sume the thin-shell approximation and ignore boundary ef¬ 
fects, the effects of mask, and the redshift dependence of the 
mean galaxy number density. While the primary effects of 
these assumptions can be corrected when making clustering 
measurements, it is likely that they couple with the geomet¬ 
ric biases considered here. Indeed, these approximations are 
likely to change the corrections due to LOS-assumptions at 
a comparable order to the corrections presented here for a 
simplified survey (although when the single-LOS method is 
used on a data with wide footprint the geometric effects are 
likely to be dominant). A more precise correction would re¬ 
quire a detailed study of how the mask effects couple with 
the geometric biases. 
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The expression in Eq. (201 is invariant under the rota¬ 


tion of coordinate system. This symmetry can be used to 
align the direction of z axis with r+ as we integrate over 
dr+. Ignoring the boundary effects and assuming that the 
redshift shell is thin, this results in in 


47r 

(2P -f 1) 


= i 

^dr-Ye,o{f-) V 


(A4) 




as the integral over d^r+ is equal to volume and only terms 
with m' = 0 because r+ is pointing in z direction. This 
expression is again invariant with respect to rotations in 
azimuthal angle (ignoring boundary effects) and can be re¬ 
duced to 


= 27ri 


47r 


(2P -b 1) 


(A5) 


d/r-Yr 0 (M-)y«o (A'-)yi!0 (mi ), 


where the m' ^ 0 terms are killed by the azimuthal in¬ 
tegral and /i- and fii are the cosines of respective po¬ 
lar angles. Since ri = r+ — r-/2 we have (1 — 

rj/i- /2)/yr — ? 7 /r + 77^/4, which will result in Eq. |22|. 


APPENDIX A: BIAS IN MOVING-LOS 
METHOD 


We will rewrite Eq. using the plane wave expansion, 
gifcr- ^ '^i‘4T^j^(^kr)Yem{k)Ye^{f), (Al) 

im 

and the addition theorem for Legendre polynomials, 

E YUk)YUf)- (A2) 

m 

This results in 

i'm' i"m" m 
Ajr 

Cr'(r-)^^W'w(r-)Y,t^,(R) 

i'^ ‘i'Kji»{kr)Ye„^„(k)Yiiirn"if-) (A3) 

A TT 


Integrating over k, by virtue of orthogonality of spherical 
harmonics, results in Eqs. (18|-(20l. 


© 0000 RAS, MNRAS 000, 
















